Role of the vortex-core energy on the Beresinkii-Kosterlitz-Thouless 

transition in thin films of NbN 
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We analyze the occurrence of the Beresinkii-Kosterlitz-Thouless transition in thin films of NbN at 
various film thickness, by probing the effect of vortex fluctuations on the temperature dependence of 
the superfluid density below Tbkt and of the resistivity above Tbkt ■ By direct comparison between 
the experimental data and the theory we show the crucial role played by the vortex-core energy in 
determining the characteristic signatures of the BKT physics, and we estimate its dependence on the 
disorder level. Our work provides a paradigmatic example of BKT physics in a quasi-two-dimensional 
superconductor. 
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Ever since the pioneering work of Berezinskii, Koster- 
litz and Thouless (BKT) 1 - 2 predicting the occurrence of 
a phase transition without a continuously broken symme- 
try in quasi 2 dimensional (2D) systems, a lot of effort has 
been devoted to study its realization in real materials 3 . 
Of particular interest have been 2D superconductors 6-13 , 
where the superconducting transition is expected to be- 
long to the BKT universality class. In these systems, 
the BKT transition can be studied through two differ- 
ent schemes. When approaching the transition tem- 
perature Tbkt from below, the superfluid density (n s ) 
(which is related to the magnetic penetration depth A) 
is expected to go to zero discontinuously at the transi- 
tion, with an "universal" relation between u s (Tbkt) and 
Tbkt itself 3 ' . Approaching the transition from above, 
one can identify the BKT transition from superconduct- 
ing fluctuations, which leave their signature in the tem- 
perature dependence of various quantities, such as resis- 
tivity or magnetization 14 . In this second scheme, the 
information on the BKT transition is encoded in the 
correlation length £(T), which diverges exponentially at 
Tbkt, in contrast to the power-law dependence expected 
within Ginzburg-Landau theory 5 . 

Many of the experimental investigations on 2D su- 
perconductors have relied on this second approach 9-12 
to identify the BKT transition through the temperature 
dependence of resistivity p(T), using eventually the in- 
terpolating formula proposed long ago by Halperin and 
Nelson 14 to describe the crossover from BKT to ordi- 
nary GL superconducting fluctuations. However, real 
superconductors have additional complicacies that can 
make such an analysis more involved than what has been 
discussed in the original theoretical approach. First, 
real systems always have some degree of inhomogene- 
ity, which tends to smear the sharp signatures of BKT 
transition compared to the clean case. As it has been 
recently shown through scanning tunneling spectroscopy 
measurements 15 ' 16 , even when disorder in the system is 



homogeneous, the system shows intrinsic tendency to- 
wards the formation of spatial inhomogeneity in the su- 
perconducting state, which has to be taken into account 
while analyzing the BKT transition. At a more funda- 
mental level, it has recently been proved experimentally 8 
that for a real superconductor the vortex core energy 
(p) can be very different from the value predicted within 
the 2D XY model, that was originally investigated by 
Kosterlitz and Thouless as the paradigmatic case to study 
the BKT transition 2 . This can give rise to a some- 
how different manifestation of the vortex physics, even 
without the extreme of a change of the order of the 
transition, as it has been proposed in the past 3,17 . Re- 
cently the relevance of p for the BKT transition has at- 
tracted a renewed theoretical interest in different con- 
text, ranging from the case of layered high-temperature 
superconductors 18-20 to the one of superconducting in- 
terfaces in artificial heterostructures 21 . 

All the above issues explain why more than 30 years af- 
ter the prediction of the BKT transition in ultrathin films 
of superconductors its occurrence in real materials is still 
controversial. The present work aims to give a paradig- 
matic example of the emergence of the BKT transition in 
thin films of NbN as the film thickness decreases. By a 
systematic comparison between A(T) and p(T), we show 
that to fully capture the 'conventional' BKT behavior in 
a real system one must account for the correct value of p 
as compared to the energy scale given by the superfluid 
stiffness J s . The analysis carried out for films of different 
thickness provides us also with an indirect measurement 
of the dependence of the vortex-core energy on disorder, 
showing that vortices become energetically more expen- 
sive as disorder increases. Such a result can be related 
to the separation between the energy scales connected to 
the SC gap (A) and J s as disorder increases, as we show 
by computing the ratio p/J s within the Bogoliubov-de- 
Gennes (BdG) solution of the attractive Hubbard model 
with on-site disorder. Our results shed new light on the 
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FIG. 1: (Color online) (a) Temperature dependence of A _2 (T) 
and p(T) for four NbN films with different thickness. The 
(black) solid lines and (red) dashed lines correspond to the 
BCS and BKT fits of the A" 2 (T) data respectively, (b) An 
expanded view of X~ 2 (T) close to Tbkt', the intersection of 
the BCS curve with the dotted line is where the BKT jump 
would be expected within the XY model, when /j, is large, 
(c) Temperature variation of R/Rn- The (red) dashed lines 
show the theoretical fits to the data, as described in the text. 



occurrence of the BKT transition on disordered films. 

Our samples consist of epitaxial NbN films grown on 
single crystalline (100) oriented MgO substrates with 
thickness (d) varying between 3-50 nm. The deposition 
conditions were optimized to obtain the highest T c ( 16K) 
for a 50 nm thick film. Details o sample preparation and 
characterisation have been reported elsewhere 8,22 . The 
absolute value of A as a function of temperature was mea- 
sured using a low-frequency (60 kHz) two-coil mutual in- 
ductance technique 8 on 8 mm diameter films patterned 
using a shadow mask. p(T) was measured on the same 
films through conventional 4-probe technique after pat- 
tering the films into 1 mmx6 mm stripline using Ar-ion 
milling. 

The first clear signature of the presence of vortices in 
our samples is provided by the deviations of A _2 (T) from 
the BCS temperature dependence as d decreases. In par- 
ticular, we observe a sharp downturn of A~ 2 (T), remi- 
niscent of the so-called universal jump of the superfluid 
density 4 . To clarify the notation, we recall that for a 2D 
superconductor J s is defined as: 
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where n 2d is the effective 2D superfluid density. In a 
conventional 3D superconductor J S (T) goes to zero con- 
tinously at the SC temperature T c . Instead within BKT 
theory the SC transition is controlled by the vortex- 
antivortex proliferation, that becomes entropically fa- 
vorable at the temperature scale Tbkt defined self- 
consistently by the relation 



ttJ s (T B kt) 



= 2. 



(2) 



Tbkt 

In the above relation the temperature dependence of 



J S (T) is due not only to the existence of quasiparticle 
excitations above the gap, but also to the presence of 
bound vortex- antivortex pairs below Tbkt- The latter 
effect is usually negligible when /i is large, as it is the 
case for superfluid films 23 . In this case one can safely 
estimate Tbkt as the temperature where the line 2T/7T 
intersects the jf CS (T) obtained by a BCS fit of the su- 
perfluid stiffness at lower temperatures. However, as /x 
decreases the renormalization of J s due to bound vor- 
tex pairs increases, and consequently Tbkt is further 
reduced with respect to T c . 8,18 To account for this effect 
we fitted the temperature dependence of A~ 2 (T) by inte- 
grating numerically the renormalization-group equations 
of the BKT theory 3 ' 18 using as only free parameter 8 the 
ratio n/J s . As input parameter for J S {T) we used the 
one obtained by a BCS fit of the data (solid lines in Fig. 
la) at low temperatures, where vortex excitations are 
suppressed, which extrapolates to zero at the mean-field 
transition temperature T c . As one can see in Fig. lb, the 
transition is still slightly rounded near Tbkt, so that the 
sharp jump is replaced by a rapid downturn at the inter- 
section with the universal 2T/tt line. We attribute this 
effect to the spatial inhomogeneity of the sample, that 
can be accounted for by assuming a distribution of lo- 
cal J](T) values around the BCS one, and performing an 
average of the A _2 (T) associated to each patch. For sim- 
plicity we assume that the occurrence probability Wi of 
each local Jj value is Gaussian, with relative width 5. We 
then rescale proportionally the local T % c and we calculate 
the resulting T BKT from the RG equations 8,19 As shown 
in Fig. la-b, such a procedure leads to an excellent fit of 
the experimental data in the whole temperature range. 
The obtained values of the ratio \ij J s (Table I) are rela- 
tively small as compared to the standard expectation of 
the XY model 24 , where 
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We recall that within the BKT approach to the XY 
model the value of /i is fixed by the cut-off at short length 
scale of the energy of a vortex line, 



E = irJ s 



log 



(4) 



where L is the system size, £o is the coherence length and 
fi = irJ s a. By mapping the (lattice) XY model into the 
continuum Coulomb- gas problem 24 one obtains a ~ 7r/2, 
so that [i attains the value in eqn. (3). However, in our 
samples p is better estimated from the loss of condensa- 
tion energy within the vortex core (see discussion below) , 
leading to a smaller p/ J s ratio and to the deviations of 
the data from the BCS fit already before the renormal- 
ized stiffness reaches the universal value 2T/n. 

To further establish the validity of the values of p ob- 
tained from the behavior of A~ 2 (T) below Tbkt we now 
use the same set of parameters to analyze the p(T) above 
Tbkt- In 2D, the contribution of SC fluctuations to the 
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conductivity can be encoded in the temperature depen- 
dence of the SC correlation length, 8a oc £ 2 (T). The 
functional form of £(T) depends on the character of the 
SC fluctuations, being power-law for Gaussian Ginzburg- 
Landau (GL) fluctuations 5 and exponential for BKT-like 
vortex fluctuations 2 ' 14 . Due to the proximity between 
Tbkt and T c (Table I) , we expect that most of the fluc- 
tuation regime for the paraconductivity will be described 
by standard GL SC fluctuations, while vortex fluctua- 
tions will be relevant only between T c and Tbkt- To 
interpolate between the two regimes we resort then to 
the Halperin-Nelson formula for £ 14 
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where t = (T — Tbkt) /Tbkt and A is a constant of 
order one. b is the most relevant parameter to deter- 
mine the shape of the resistivity above the transition 
and is connected 21 both to the relative distance t c be- 
tween Tbkt and T c , t c = (T c ~ T B kt)/T B kt, and to 
the value of p: 
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The normalized resistance corresponding to the SC cor- 
relation length (5) is given by 
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where Rn is the normal-state resistance (that we take 
here as Rn = R(T — 1.5Tbkt)- Finally, to account 
for sample inhomogeneity, we map the spatial inhomo- 
geneity of the sample in a random-resistor-network prob- 
lem, by associating to each patch of stiffness J* a nor- 
malised resistance p.; = Ri/ Rn obtained from Eq. (7) by 
using the corresponding local values of T l c and T % BKT com- 
puted above. The overall sample normalised resistance 
p = R/Rn is then calculated in the so-called effective- 
medium-theory (EMT) approximation 25 , where p is the 
solution of the self-consistent equation 
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and Wi is the occurrence probability of each resistor, i.e. 
of the corresponding J* value, as determined by the anal- 
ysis below Tbkt- As it has been discussed in Ref. 26 , the 
EMT approach turns out to be in excellent agreement 
with the exact numerical results for a network of resis- 
tors undergoing a metal-superconductor transition, even 
in the presence of SC fluctuations. We can then em- 
ploy Eq. (8) to compute R/ Rn of our samples, by using 
the probability distribution of width S known from the 
analysis of A(T), and by treating A and b as free param- 
eters. The resulting fits are in excellent agreement with 
the experimental data (Fig. lc). Moreover, considering 
that the interpolation formula (5) between the BKT and 



TABLE I: Sheet resistance (R s ), Magnetic penetration depth 
(A(T —¥ 0)), Tbkt, BCS transition temperature T c , along 
with the best fit parameters (see text) obtained from the BKT 
fits of A~ 2 (T) below Tbkt and of R(T) above Tbkt for NbN 
thin films of different thickness d. The temperature T c is 
obtained by the extrapolation of the BCS fit of A -2 well below 
Tbkt- 
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FIG. 2: (Color online) (a) Numerical results for the disorder 
dependence of p/J s and A/J s as a function of disorder for 
the attractive Hubbard model, (b) Experimental values for 
the same ratios in our NbN films, plotted as a function of the 
normal-state sheet resistance Rs- 



GL fluctuation regime is necessarily an approximation, 
the obtained values of b are in very good agreement with 
the theoretical estimate (6) (Table I). Thus, our analysis 
above T c not only provides us with a remarkable exam- 
ple of interpolation between the GL and BKT fluctuation 
regimes, but it also demonstrates the validity of the val- 
ues of p obtained from A(T). Finally, we would like to 
stress that b cannot be used completely as a free param- 
eter while fitting the p(T). Attempting to fit the BKT 
fluctuation regime at T^>Tbkt (as proposed in the liter- 
ature [10,11]), results in unphysical b values with respect 
to relation (6). 

Once established the robustness of our estimate of p, 
we discuss now the values reported in Table I, and their 
thickness dependence. We first notice that the values of 
p obtained by our fit are of the order of magnitude of the 
standard expectation for a BCS superconductor. Indeed, 
in this case one usually 27 estimates the p as the loss in 
condensation energy within a vortex core of size of the 
order of the coherence length £o> 
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is the condensation-energy density. In the 
clean case Eq. (9) can be expressed in terms of J s by 
means of the BCS relations for e cone i and £q. Indeed, since 
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tcond = iV(0)A 2 /2, where N(0) is the density of states at 
the Fermi level and A is the BCS gap, and £o = £bcs — 
Hvf/ttA, where vf is the Fermi velocity, and assuming 
that n s = n at T = 0, where n = 2N(0)v F m/3, one has 
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so that it is quite smaller than in the A"Y"-model case 
(3). While the exact determination of n depends on small 
numerical factors that can slightly affect the above esti- 
mate, the main ingredient that we should still account 
for is the effect of disorder, that can alter the relation 
between e con( j, A and J s and explain the variations ob- 
served experimentally. To properly account for it we 
computed explicitly both fi and J s within the attractive 
two-dimensional Hubbard model with local disorder: 

H = -t^2 C L C j" + h - c - - \ U \ X] n ^ n ^ + X! V i n iai 
(ij)a i i<7 

(11) 

which we solve in mean field using the BdG equations 28 . 
The first sum is over nearest-neighbors pairs and we work 
on a N = N x x N y system, with a local potential Vi ran- 
domly distributed between < Vi < Vq. J s is computed 
by the change in the ground-state energy in the presence 
of a constant vector potential 29 , while /i is computed by 
means of Eq. (9), by determing both e conc i and £ in the 
presence of disorder 30 at doping n = 0.87 and coupling 
U/t = 1. The resulting value of fi/J s at T = is reported 
in Fig. 2a: It is of order of the BCS estimate and it shows 
a steady increase as disorder increases, in agreement with 
the experimental results, shown in Fig. 2b, where we take 
the normal state sheet resistance R K as a measure of dis- 



order as the film thickness decreases. This behavior can 
be understood as a consequence of the increasing separa- 
tion with disorder between the energy scales associated 
respectively to the A, which controls e con d, and J s , as 
it is shown by the ratio A/J s that we report in the two 
panels of Fig. 2 for comparison. Notice that even though 
we used a weaker coupling U/t = 1 as compared to other 
recent studies 29,31 this is still a large coupling strength 
as compared to our NbN samples, so that the numeri- 
cal values of A/J s are larger than experimental ones 32 . 
Nonetheless, our approach already captures the experi- 
mental trend of \xjJ s as a function of disorder, and its 
correlation with the A/J s behavior at large disorder. 

In summary, we have shown that to correctly identify 
the typical signatures of the BKT transition in thin films 
of NbN we must properly account for /i values smaller 
than expected within the standard approach based on 
the XY model 1 ' 2 . We also observe a steady increases of 
the ratio n/J s as the film thickness decreases. This effect 
can be understood within a model for disordered super- 
conductors, resulting from increasing separation between 
the energy scales associated with A and J s . It would be 
interesting to investigate if a similar effect could be at 
play also in other systems, as disordered films of InO^ 7 , 
or high-temperature cuprate superconductors, where a 
large n value has been indirectly suggested by the anal- 
ysis of the superfluid density data 18 . 
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